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Abstract: 



o ■ 

Q\ ■ It is shown that the nonlinear Ermakov-Milne-Pinney equation p + v(x)p = a/p obeys the 
property of covariance under a class of transformations of its coefficient function. This property 
is derived by using supersymmetric, or Darboux, transformations. The general solution of the 
transformed equation is expressed in terms of the solution of the original one. Both iterations of 
these transformations and irreducible transformations of second order in derivatives are consid- 
ered to obtain the chain of mutually related Ermakov-Milne-Pinney equations. The behaviour 
of the Lewis invariant and the quantum number function for bound states is investigated. This 
construction is illustrated by the simple example of an infinite square well. 
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1. The interrelation between a special nonlinear differential equation, known as the Er- 
makov, Milne or Pinney [ |lj equation (EMP equation), and a corresponding linear equation 
has been investigated in great detail by many authors (see [ ||, the recent review in [ [| and 
references therein). 

In classical dynamics these equations appear most often in context with time-dependent os- 
cillators in the form 

p + uj 2 (t)p = ap~ 3 (1) 
and 

q + LU 2 (t)q = 0, (2) 
where the so-called Lewis invariant [ HI 
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a ^2 + O 9 ? ~ fa) ( 3 ) 
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plays a role as a constant of motion for (|2j). For a recent application to pulse induced transitions 
see[§. 

In Quantum Mechanics, the general solution of the Schrodinger equation 

1 p" + k\x)^ = , k\x) = ^{E-V{x)) (4) 
can be expressed in amplitude-phase form as 

i/j(x) = a p(x) sin ( / p~ 2 (x') dx' + /3j , (5) 
where a, f5 are constants and p(x) is an arbitrary solution of 

p" + k 2 (x)p = ap- 3 . (6) 
Eq. (|5|) directly implies the condition 

i r +ao i 

N(E) = - p- 2 (x,E)dx = n + 1 , n= 0,1,2,... (7) 

7T V-oo 



for the bound state energies E n . 

In reverse, the Ermakov method can be considered as a method to solve the EMP equation. 
Here we will discuss some aspects of the Darboux [ 0], or super symmetric [ |7|], covariance of 
this nonlinear differential equation. The role of Darboux transformations in the study of soliton 
solutions of nonlinear evolutionary equations, like KdV and KP, is well known (see, for example, 
[ P| and references therein). Nevertheless, up to our knowledge, this approach was not used for 
an investigation of EMP equation. 

2. We start from the standard one-dimensional Schrodinger equation (H = l,m = 1/2): 
- i//'(x, E) + (W 2 (x) - W'(x))il>(x, E) = Ei/>(x, E), (8) 



where for convenience the potential V(x) is written in "supersymmetric form" [ [7j in terms of 
the real superpotential W(x) : 

V(x) = W 2 (x) - W'{x). (9) 

The Hamiltonian H = —d 2 + V(x) is intertwined with another (superpartner) Hamiltonian 

H = -d 2 + V{x); (10) 

V{x) = W 2 {x) + W'{x) (11) 

by the component of supercharge 

q + = -d + W(x); (12) 

Hq + = q + H. (13) 
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This intertwining relation, together with the hermitian conjugated one, 

q-R = Hq- (14) 
q- = ( 9 +)t = d + W(x), (15) 

seem to be the most important elements of standard super symmetrical Quantum Mechanics 
(SUSY QM) [ |?J and its generalizations [ [9|, [l(], [ll], 12]. In particular, just these relations ([H 



( |j~4|) lead to the isospectrality (possibly, up to zero modes of operators q^) of Hamiltonians H, H 
and to the connection between their normalized eigenfunctions, if E belongs to a discrete part 
of the spectrum: 

i>(x,E) = ^=q-i>(x,E) (16) 
V E 



^(x,E) = ^=q + ^(x,E). (17) 



Thus, if the model with Hamiltonian H is exactly solvable (or quasi-exactly-solvable [ [l3|]), the 
model with Hamiltonian H is also exactly solvable (or quasi-exactly-solvable). By iteration of 

this procedure, from V = W 2 + W = W 2 — W 1 + c; c = const to the next potential V = 
W 2 + W' + c, one can obtain a chain of isospectral (up to zero modes of q^) models with 
known eigenvalues and eigenfunctionsQ In terminology of Mathematical Physics we described 
now schematically the well known Darboux-Crum [ |l4|] transformations [ § for Sturm-Liouville 
operator. 

Let us denote two pairs of linearly independent solutions (may be both in the pair are 
not normalizable) of the stationary Schrodinger equations with Hamiltonians H and H by 
ipi(x,E), ip2(%,E) and ipi(x,E), ip2(%,E), correspondingly. By using ([i~6|), (|i~7D, one can check 
that their Wronskians coincide: 

A = if>iip2 — ^1^2 = Wi^2 — ipiV2 = A = const. (18) 



As was found [ 15], an arbitrary pair of linearly independent solutions ipi, ip2 of linear equation 



(B) leads to a general solution 



p(x,E) 



1/2 

Aipf (x, E) + B^ 2 {x, E) + 2Ctyi(s, E)ifo(x, E) 



(19) 



of the nonlinear EMP equation (for simplicity we normalized p in Eq.(|]) as a = 1): 
- p"(x,E) + (w 2 (x) - W'(x,E) - e) P (x,E) = (20) 
where A, B, C above are arbitrary constants, which have to satisfy the relation 
AB-C 2 = ±. (21) 

1 In notations of [ Eel this is a class of so called reducible second order supertransformations with a shift c. 
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By a similar procedure with an arbitrary pairQ i/Ji(x, E), ^{x, E) of solutions of the partner 
Schrodinger equation with Hamiltonian (|10|), ([ll]), one can construct a general solution 



p(x,E) 



A${x,E) + B$ 2 {x,E) + 2C$ 1 (x,E)ifo(x,E) 



1/2 



(22) 



of the partner EMP equation (with the partner coefficient function (V(x) — E)): 
- p"(x, E) + ( W 2 (x) + W'{x, E) - E] p(x, E) 



p{x,E) 



(23) 



After substitution of (|l6|) in (p2|), one obtains a general solution p(x,E) of (p3|), which is 
expressed now not in terms of ip\, ip2, but in terms of linearly independent solutions ip\, tp2 of 
the partner Schrodinger equation (j8|). In such a way we found an indirect connection between 
solutions p(x, E) and p(x, E) of two EMP equations via direct connection of solutions of "their" 
Schrodinger equations. But this prognosticated construction would become much more elegant, 
if the general solution p(x, E) of ( p3| ) could be expressed directly in terms of solution p(x, E) 
of (p0|). Indeed, this is the case: straightforward calculations lead to the following nonlinear 
expression: 



Ep 2 (x,E) = (p'(x,E)) +W(x)^p 2 (x,E^j +W 2 (x)^p(x,E)^j + — 1 



E) 



d + W{x) )p(x,E) 



1 2 



+ 



p 2 (x,EY 



(24) 



So, the nonlinear EMP equation possesses the following property, similar to the supersymme- 
try of Schrodinger equation (or to Darboux covariance of Sturm-Liouville operator): if p(x,E) 



satisfies equation (20), its nonlinear transformation (p4|) satisfies the same EMP equation, how- 
ever with a coefficient function transformed accordingly to (|23|). 

By iterations one obtains again a chain of such EMP equations with related coefficient func- 
tions ("potentials"), whose general solutions can be expressed consecutively in terms of the so- 
lution p(x, E) of the first original equation (|20|). This is a nonlinear analogue of Darboux-Crum 
reducible transformations for the Sturm-Liouville equation. In particular, for two consecutive 
transformations: 



V = W 2 - W' -» V = W 2 + W' = W 2 - W' + c -c V = W 2 + W' + c: 



p^ p^ p, 



(25) 



the solution p of EMP equation ( p0[ ) with coefficient function V can be expressed in terms of 
original solution p : 



E 2 f 



w(w + W S - E 
(W + W) 2 



2 / \ 2 

.2 i / w i W 1 / J\2 



p 2 + [w + w) {ply + iw + w 



w\w + w S - L 



(p 2 )' 

(26) 



We choose here the same constants A, B, C and take into account the equality (N M between Wronskians 
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where W(x) and W(x) are connected by equality (|2 
single function f(x): 



and can be expressed [ 1C] in terms of a 



W = f 



2f 



4/ 



W = f + 



2f> 



The transformation 



4/ 

can be inverted: 



Ep 2 (x,E) 



d-W(x) )p(x,E) 



+ - 



1 



p{x,E\ 



Though the difference in signs in front of W(x) in (12 
relation between (||) and QlT|), the direct check of 
the intermediate auxiliary equality for p and p : 

Wp 2 . 



(27) 



(28) 

and (|28|) seems to be obvious due to the 
) is not so trivial: it is useful to explore 

(29) 



pp' - Wp = -pp' 

3. The various generalizations of standard SUSY QM transformations were investigated in the 
literature. In particular, SUSY transformations of higher order in derivatives were constructed 
in 



1C, 11 1 (see also the recent papers [ 16]). It was shown 



10 1 for the one-dimensional 



Schrodinger equation, that only two "elementary" types of transformations exist: the first order 
ones, generated by q^ ( see (12), ( |i~5| ) ), and the irreducible second order transformations, 
which can not be expressed as a product of two standard first order transformations. All 
supercharges (and supertransformations) of higher order can be represented as a number of first- 
and second-order steps [ 10|. In the context of the present paper it seems to be interesting, in 
addition to the construction which lead to ( |24|) and (28), to take into account these second order 
supertransformations of Schrodinger operator and their possible analogues for EMP equation. 

The most general form of solution of intertwining relations (|l~3| ) , but with second order super- 
charges , was found in [ [T(| : 

q + = d 2 -2f(x)d + b(x); (30) 
q- = ( 9 +)t =d 2 + 2f{x)d + 2f{x)+b{x)- (31) 



b(x) 
V(x) 
V(x) 



-2f'(x 
+2f'(x 



+ f(x) + 
+ f\x) + 



f"(x) , f fjx) 
2/0*0 



Vix) + 
f"(x) 



+ 



d 



4/2 (x) 



Vix) 
fix) 



2/(x) 
fix) 



d 



2f(x) \2f(x)J 4/ 2 (x)' 



4/2 ( X ) 
d 



(32) 
(33) 
(34) 



where f(x) is an arbitrary real function. For reducible transformations d = —c 2 /4, it was 



introduced already in Eq.(|27|), and an irreducible situation corresponds to the case when the 
constant d is positive. 
Substitution ofQ 

$i{x,E) = -=±==q-4i(x,E)l » = 1,2 (35) 
VE Z + a 



3 The normalization factor is obtained from the corresponding relation of the algebra of Second Order SUSY 
Quantum Mechanics [ §, 0]: q + q~ = H 2 + d; d > 0. 
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with the second order q~ from (|31|) , (|32|), into Eq.(|22|) leads to new expression for p(x, E) in terms 
of original solutions ipi(x,E), ^2(x,E). This expression, after replacing ^"s via Schrodinger 
equation @, again can be written in terms of function p(x, E) and its first derivative only: 

(E 2 + d)p 2 (x,E) = Af(x)L'(x,E)\ + ^2f 2 (x)-f'(x)-E^p 2 (x,E) + 



+ 2f(x)(2f 2 (x)-f'(x)-E^p 2 (x,E) 



4. One of the reasons of the rather long peculiar interest to Ermakov systems and to related 
EMP equations is connected with existence of the x— independent quantity (or t— independent in 
the context of classical harmonic oscillator instead of Schrodinger equation (||)). This invariant, 
so called Ermakov-Lewis invariant, is expressed in terms of an arbitrary solution fy(x, E) of (|8|) 
and an arbitrary solution p(x,E) of the related "auxiliary" equation (20): 

i>(x,Ey " 



7 = 1 
2 



p(x,E) 



+ p(x, E)ip'(x, E) - p'(x, E)ij}(x, E) 



(37) 



Calculation of ( |37| ) with p and ^ instead of p and ^ shows that the value of Ermakov-Lewis 
invariant is not changed under this supersymmetric transformation and its iterations, i.e. 1 = 1. 

The quantum number function (|7|) for bound states E n transforms to an analogous function 
N(E n ) for the system (|To|). N(E n ) either coincides with N(E n ) or differs by ±1. This difference 
depends on the interrelation of spectra H and H, i.e. [ [l^] on the asymptotic properties of 
superpotential W(x). For the case of second order transformations (|3l| ) this difference can be 
equal to ±2 (see [0]). 

5. As an illustration, we will discuss a simple example in some detail, namely the infinite 
square well potential 

f -1 \x\ < vr/2 

V(x) = ' , (38) 

I co \x\ > vr/2 

which can be expressed (^) in terms of the superpotential 

[ tanx Ixl < 7r/2 

W(x) = (39) 

[ OO \X\ > 7T/2 

In the following we will confine ourselves to the well region |x| < ir/2 unless otherwise stated. 
The supersymmetric partner is: 

V(x) = W 2 (x) + W'{x) = -l + 2sec 2 x (40) 
for the supersymmetric partner potential. The eigenvalues for V{x) are 

E n = n{n + 2), n = 0,1,2,..., (41) 

where by construction the ground state eigenvalue Eq is zero. Because of super symmetry, the 
eigenvalues coincide with those of V{x): 

E n+1 = E n , n = 0,l,2,.... (42) 
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According to Eq. fll9|) we can now write the general solution of the EMP-equation (^Q) 

p"(x) + k 2 p(x) = p- 3 {x), k 2 = E + l (43) 

in terms of the solutions ipi(x) = sin kx and ip2(%) = coskx of the Schrodinger equation: 

p 2 (x) = A sin 2 kx + B cos 2 kx + 2C sin kx cos kx , AB - C 2 = l/k 2 . (44) 

Note that p{x) is infinite outside the potential well. Assuming the special case C = 0, A = B = 
l/k for simplicity, we have the solution 

p 2 (x) = l/k = const. (45) 
As a check of the Milne quantization condition (|7|) one obtains 

N(E) = - p- 2 dx = - kdx = k = n + l, n = 0,1,2,... (46) 

J-oo K J-n/1 

which reproduces, of course, Eq.(|4l|) 

By means of the super symmetric covariance we can now immediately write down the solution 
(H) of the EMP-equation (H) for the partner potential @: 

p 2 (x) = (p' + ptanx) 2 + p~ 2 = k~ 1 tm 2 x + k = ~f(k 2 + tan 2 x) (47) 
with 7 = 1/kE. 

By some elementary algebra, one can directly check that p is indeed a solution of the EMP 
equation (23). 

From the solution ([17]) we can also compute the quantum number function 
1 f+°° Ar 1 7r 

^-Lr^^- 1 <48) 



and with N(E) = k from ( jig ) we arrive at the relation N(E) = N(E) + 1 as was predicted, 
taking into account Eq. (p2]) . 

Finally, we can also verify that the values of the invariants I and I agree. If ^f(x) is any 
solution of the Schrodinger equation for the square well potential V{x) and p 2 (x) = l/k is taken 
from (fH|), we have 

21 = + (p^' - p'mf = jfe*g + k- 1 ^' 2 . (49) 

The partner wave function is \[E ^ = + tan x ^ with \[E = (tan 2 x — E) ^ + tan x \E' / . We 
evaluate the (constant) value of the invariant / at x = 0, which gives with ^(0) = etc- an d 

21 = (^V^^-p^) 2 
v Po ' 

^ 2 1 ^- \2_ t -l lTf /2 , i >lT ,2 



7^ 2 



+ 7 A; 2 (- v^o) =fc-^^ + A;^ = 2/. (50) 
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